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Considering gravitino dark matter scenarios, we obtain tlie full gauge-invariant result for the relic 
density of thermally produced gravitinos to leading order in the Standard Model gauge couplings. 
For the temperatures required by thermal leptogenesis, we find gaugino mass bounds which will be 
probed at future colliders. We show that a conceivable determination of the gravitino mass will 
allow for a unique test of the viability of thermal leptogenesis in the laboratory. 
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INTRODUCTION 

The smallness of the neutrino masses can be under- 
stood naturally in terms of the see-saw mechanism [l| 
once the Standard Model is extended with right-handed 
neutrinos which have heavy Majorana masses and only 
Yukawa couplings. For a reheating temperature after in- 
flation, Tr , which is larger or not much smaller than the 
masses of the heavy neutrinos, these particles are pro- 
duced in thermal reactions in the early Universe. The 
CP-violating out-of-equilibrium decays of the heavy neu- 
trinos generate a lepton asymmetry that is converted into 
a baryon asymmetry by sphaleron processes This 
mechanism, known as thermal leptogenesis, can explain 
the cosmic baryon asymmetry for Tr > 3 x 10^ GeV 3]. 

One will face severe cosmological constraints on Tr if 
supersymmetry (SUSY) is discovered. An unavoidable 
implication of SUSY theories including gravity is the ex- 
istence of the gravitino G which is the gauge field of local 
SUSY transformations. As the spin-3/2 superpartner of 
the graviton, the gravitino is an extremely weakly inter- 
acting particle with couplings suppressed by inverse pow- 
ers of the (reduced) Planck scale Mr = 2.4 x lO^^ GeV [1]. 
In the course of spontaneous SUSY breaking, the grav- 
itino acquires a mass and the couplings of the spin- 
1/2 goldstino which become dominant for small mg. De- 
pending on the SUSY breaking scheme, mg can range 
from the eV scale up to scales beyond the TeV region [5j . 
Gravitinos can be produced efficiently in the hot primor- 
dial plasma. Because of their extremely weak interac- 
tions, unstable gravitinos with < 5 TeV have long 
lifetimes, > 100 sec, and decay after big bang nucleo- 
synthesis (]3BN). The associated decay products affect 
the abundances of the primordial light elements. De- 
manding that the successful BBN predictions are pre- 
served, bounds on the abundance of gravitinos before 
their decay can be derived which imply Tr < 10^ GeV 
for TTig < 5 TeV fij. Thus, the temperatures needed for 
thermal leptogenesis are excluded. 

We therefore consider SUSY scenarios in which a grav- 
itino with ttIq Z 10 GeV is the lightest supersymmetric 
particle (LSP) and stable due to R-parity conservation. 



These scenarios are particularly attractive for two rea- 
sons: (i) the gravitino LSP is a compelling dark matter 
candidate and (ii) thermal leptogenesis can still be a vi- 
able explanation of the baryon asymmetry [7|. 



THERMAL GRAVITINO PRODUCTION 

Assuming that inflation governed the earliest moments 
of the Universe, any initial population of gravitinos must 
be diluted away by the exponential expansion during the 
slow-roll phase. We consider the thermal production (or 
regeneration) of gravitinos in the radiation-dominated 
epoch that starts after completion of reheating at the 
temperature Tr. Gravitinos with 10 GeV are not 

in thermal equilibrium with the primordial plasma after 
inflation because of their extremely weak interactions. At 
high temperatures, gravitinos are generated in scattering 
processes of particles that are in thermal equilibrium with 
the hot SUSY plasma. The calculation of the relic den- 
sity of these thermally produced gravitinos, , requires 
a consistent finite-temperature approach. A result that 
is independent of arbitrary cutoffs has been derived for 
SUSY quantum chromodynamics (QCD) in a gauge in- 
variant way [8]. Following this approach, we provide the 
complete SU(3)c xSU(2)l xU(1)y result to leading order 
in the couplings. 

We compute ilj^ from the Boltzmann equation for the 



gravitino number density 



dt 



3Hn^ — Cq 



(1) 



The term proportional to the Hubble parameter H ac- 
counts for the dilution by the cosmic expansion. The 
collision term Cq describes the production and disap- 
pearance of gravitinos in thermal reactions in the pri- 
mordial plasma. Since the phase space density of the 
gravitino is significantly below those of the particles in 
thermal equilibrium, gravitino disappearance processes 
can be neglected. Thus, Cg is given by integrating the 
thermal gravitino production rate. 

Considering a primordial plasma with the particle con- 
tent of the minimal SUSY Standard Model (MSSM) in 



2 



the high-temperature hmit, we calculate the thermal 
production rate of gravitinos with E > T using the 
Braaten-Yuan prescription Q and hard thermal loop 
(HTL) resummation [lO| . With the systematic treatment 
of screening effects in the plasma, we obtain a finite re- 
sult in a gauge-invariant way. Moreover, in contrast to 
previous estimates [12| , our result does not depend on 
arbitrary cutoffs. The explicit form of the thermal grav- 
itino production rate and the detailed derivation will be 
presented in a forthcoming publication 13 1. 

After numerical integration of the thermal gravitino 
production rate, we obtain the SU(3)cXSU(2)lxU(1)y 
result for the collision term: 
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where the gaugino mass parameters Mi, the gauge cou- 
plings gi, and the constants ci and ki are associated with 
the gauge groups U(1)y, SU(2)l, and SU(3)c as given in 
Table [H In expression ^ the temperature T provides 
the scale for the evaluation of Mi and gi. Note that HTL 
resummation [9, J_0J requires weak couplings, <C 1, and 
thus high temperatures T ;» 10^ GeV. 

Our result = 1.271 for the SU(3)c contribution is 
larger than 1.163 obtained from 3J- This results from an 
analytical disagreement [l3|: We find a cancellation of 
the term T^{N + nf)[U2{-e~^^'^) -tt^/G] given as part 
of /bfb in (C.14) of Ref. d. 

Assuming conservation of entropy per comoving vol- 
ume, the Boltzmann equation (jlj can be solved analyti- 
cally [i,[lil. With the coUision term (g]), we find 
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with the Hubble constant h in units of 100 kmMpc^ s^^ 
and the constants Ui given in Table HI Here Mi and gi 
are understood to be evaluated at Tr. With our new 
^3 value, wc find an enhancement of about 30% of the 
SU(3)c contribution to the relic density. 
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1: The relic gravitino density from thermal production, 
as a function of Tr. The solid and dashed curves show 
the SU(3)cXSU(2)lxU(1)y results for universal (Mi, 2,3 = 
mi/2) and non-universal (0.5 Mi, 2 = M3 = mi/2) gaugino 
masses at Mgut, respectively. The dotted curves show our 
result of the SU(3)c contribution for M3 = mi/2 at Mqut. 
The gray band indicates the dark matter density f2DM/i^. 



Figure [T] shows f2~^ft.^ as a function of for = 10, 
50, and 300 GeV. With mi/2 = 400 GeV, the sohd 
and dashed lines are obtained respectively for universal 
(Afi,2,3 = "11/2) and non-universal (0.5 Mi, 2 = M3 = 
TO1/2) gaugino masses at the grand unification scale 
Mgut 2 x 10^'^ GeV. The SU(3)c contributions are 
shown by the dotted lines. The gray band indicates the 
dark matter density [l5{ 



r^DM/i = 0.105 



+0.007 
-0.010 



(4) 



We find that electroweak processes enhance il- by 
about 20% for universal gaugino masses at Mqut- In 
non-universal cases. Mi. 2 > M3 at AfcuT, the electro- 
weak contributions are more important. For 0.5Afi,2 — 
M3 at Mgut, they provide about 40% of Vl^-^ . 



TABLE I: Assignments of the index i, the gauge coupling 
gi, and the gaugino mass parameter Mi to the gauge groups 
U(1)y, SU(2)l, and SU(3)c and the values of the associated 
constants Ci, ki, and Ui. 
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COLLIDER PREDICTIONS OF LEPTOGENESIS 

Thermal leptogenesis requires Tr si 3 x 10^ GeV 
This condition together with the constraint flj^ < Qdm 
leads to upper limits on the gaugino masses. The SU(3)c 
result for implies limits on the gluino mass [^, [3] . 
With our SU(3)cXSU(2)lxU(1)y result, the limits on 
the gluino mass M3 become more stringent because of 
the new k^ value and the additional electroweak contribu- 
tions. Moreover, as a prediction of thermal leptogenesis. 
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FIG. 2: Upper limits on the gaugino mass parameter mi/2 
from ti^-' < r^DM for the indicated values of Tr. The solid 
and dashed curves show our SU(3)cXSU(2)lxU(1)y results 
for universal (Mi, 2, 3 = "^1/2) and non-universal (0.5 Mi, 2 = 
M3 = mi/2) gaugino masses at Mgut, respectively. The 
dotted curves show the SU(3)c limits for A/3 — mi/2 at Mqut. 



we obtain upper limits on the electroweak gaugino mass 
parameters Mi, 2. At the Large Hadron CoUider (LHC) 
and the International Linear Collider (ILC), these limits 
will be probed in measurements of the masses of the neu- 
tralinos and charginos, which are typically lighter than 
the gluino. If the superparticle spectrum does not re- 
spect these bounds, one will be able to exclude standard 
thermal leptogenesis. 

Figure [2] shows the gaugino mass bounds for Tr — 10^, 
3 X 10^, and 10^° GeV evolved to Mqut, i-e., in terms 
of limits on the gaugino mass parameter mi/2- Here 
^DM = 0.126 is adopted as a nominal Scr upper limit 
on r^DM^^- With the observed superparticle spectrum, 
one will be able to evaluate the gaugino mass parameters 
Ml, 2, 3 at Mqut using the SUSY renormalization group 
equations 13, 3 While the determination of Mi, 2 
at low energies depends on details of the SUSY model 
that will be probed at coUiders Q, the bounds shown in 
Fig. [2] depend mainly on the Mj relation at Mqut- This 
is illustrated by the solid and dashed curves obtained 
with Ml, 2, 3 — TOi/2 and 0.5 Mi, 2 = M3 = mi/2, respec- 
tively. The dotted curves represent the SU(3)c limits for 
M3 = TOi/2 at Mqut and emphasize the importance of 
the electroweak contributions. 



DECAYS OF THE NEXT-TO-LIGHTEST 
SUPERSYMMETRIC PARTICLE 

With a gravitino LSP of > 10 GeV, the next- 
to-lightest SUSY particle (NLSP) has a long lifetime of 



tnlsp ~ 10 s [20|, |2l| . After decoupling from the primor- 
dial plasma, each NLSP decays into one gravitino LSP 
and Standard Model particles. The resulting relic density 
of these non-thermally produced gravitinos is 



"NLSP" 



'7^NLSP 



(5) 



where ttt-nlsp is the mass of the NLSP and riNLSP^^ is the 
relic density that the NLSP would have today, if it had 
not decayed. As shown below, more severe limits on TO1/2 
are obtained with VL^^h? taken into account. Moreover, 
since the NLSP decays take place after BBN, the emit- 
ted Standard Model particles can affect the abundance 
of the primordial light elements. Successful BBN pre- 
dictions thus imply bounds on and tonlsp 23, 21|. 
From these cosmological constraints it has been found 
that thermal leptogenesis remains viable only in the cases 
of a charged slepton NLSP or a sneutrino NLSP [l^, [l^l . 

Note that the cosmological constraints from BBN can 
become much weaker with entropy production after de- 
coupling of the NLSP and before BBN. For example, 
large parts of the parameter region disfavored by BBN 
constraints on charged slepton NLSP scenarios become 
allowed with a moderate amount of entropy produc- 
tion (23I I. In addition, such an entropy production dilutes 
both the generated baryon asymmetry and the thermally 
produced abundance of gravitinos. Therefore, upper lim- 
its on the gaugino masses will still allow us to probe the 
viability of thermal leptogenesis. However, the Tr labels 
of the curves in Fig. [2] will change to higher values. 



COLLIDER TESTS OF LEPTOGENESIS 

Thermal leptogenesis will predict a lower bound on 
the gravitino mass rriQ once the masses of the Standard 
Model superpartners are known. With a charged slepton 
as the lightest Standard Model superpartner, it could 
even be possible to identify the gravitino as tlie_LSP and 

11, 



24 



to measure its mass mg at future colliders 
Confronting the measured with the predicted lower 
bound will then allow us to decide about the viability of 
thermal leptogenesis. 

In order to explain our method, we do now con- 
sider an exemplary SUSY model. Let us assume that 
the analysis of the observed spectrum [l^ will point 
to the universality of the soft SUSY breaking param- 
eters at Mgut and, in particular, to the minimal su- 
pergravity (mSUGRA) scenario with the gaugino mass 
parameter mi/2 — 400 GeV, the scalar mass parameter 
Too = 150 GeV, the trilinear coupling Aq — —150 GeV, 
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FIG. 3: Probing the viability of thermal leptogenesis. The 
solid curves show the limits on the gaugino mass parameter 
mi/2 from + n^'^P < ng^" for Tr = 10^ 3 x 10^ and 
10^° GeV. The dashed line indicates the mi/2 value of the 
considered scenario. The vertical solid line is given by the ri 
NLSP mass which limits the gravitino LSP mass from above: 
mg < — 143.4 GeV. The values at which tempera- 
tures above 3 x 10^ GeV and 10*' GeV remain allowed are in- 
dicated by the dark-shaded (dark-green) and medium-shaded 
(light-green) regions, respectively. The mg values within the 
light-shaded region are excluded by BBN constraints. 



Figure [3] allows us to probe the viability of thermal 
leptogenesis in the considered mSUGRA scenario.^ From 
the constraint H-^ -I- < ri^^, we obtain the solid 

curves which provide the upper limits on mi/2 for = 
10^ 3 X 10^ and 10^° GeV. The dashed line indicates the 
r7ii/2 value of the considered scenario. The vertical sohd 
line is given by mf^ — 143.4 GeV which limits ttiq from 
above. In the considered scenario, the mi/2 value exceeds 
the limits for > 10^" GeV. Thus, temperatures 

above 10^° GeV can be excluded. Temperatures above 
3 X 10^ GeV and 10^ GeV remain allowed for rriQ values 
indicated by the dark-shaded (dark-green) and medium- 
shaded (light-green) regions, respectively. The values 
indicated by the light-shaded region are excluded by BBN 
constraints for late ri NLSP decays."^ 

Here thermal leptogenesis, Tr > 3 x 10^ GeV, pre- 
dicts ^ 130 GeV and thus a ti lifetime of Tf-^ > 
10^^ s [20|, |21[. If decays of long-lived ti's can be 



analyzed at colliders giving evidence for the gravitino 
LSP [23, m, m, there will be the possibility to deter- 
mine mg in the laboratory: From a measurement of the 
lifetime Tfj^ governed by the decay ti — > Gr, mg can be 
extracted using the supergravity prediction for the asso- 
ciated partial width. 



487rTO|M| 



(6) 



as obtained for mr — > 0. Moreover, for 
TOq can be infered kinematically from the energy of the 
tau, Er, emitted in the 2-body decay ti Gt 24, 2 



a positive higgsino mass parameter, fj, > 0, and the mix- 
ing angle tan/3 = 30 in the Higgs sector. A striking 
feature of the spectrum will then be the appearance of 
the lighter stau ti with m^j — 143.4 GeV as the light- 
est Standard Model superpartner [27] . In the considered 



gravitino LSP case, 10 GeV 



TOg < TOf, 



this stau is 



the NLSP and decays with a lifetime of Tf ^ > 10 s into 
the gravitino. For the identified mSUGRA scenario and 
the considered reheating temperatures, the cosmological 
abundance of the ri NLSP prior to decay can be com- 
puted from r^NLSP^^ = — 3.83 x 10~^, which is 
provided by the computer program micrOMEGAs [2^ . 
For given m^, this abundance determines il^^h^ and 
the release of electromagnetic (EM) and hadronic en- 
ergy in Ti NLSP decays governing the cosmological con- 
straints sm.! 



^ Considering bound-state effects of long-lived negatively charged 
particles on BBN, it has recently been claimed that the stau 



NLSP abundance for ■ 



10'^ 



- 10 s is severely constrained by 



the observed ^'^Li abundances [29|, |30||. If the associated bounds 
are confirmed, the considered mSUGRA scenario will be cosmo- 



m: 



? - 2mf, Er 



(7) 



While within the dark-shaded (dark-green) region 
will favor thermal leptogenesis, any mg outside of the 
medium-shaded (light-green) region will require either 
non-standard mechanisms lowering the Tr value needed 
for thermal leptogenesis or an alternative explanation of 
the cosmic baryon asymmetry. 



logically allowed only with late-time entropy production |23l |. 
Thermal leptogenesis requires right-handed neutrinos and thus 
an extended mSUGRA scenario. This could manifest itself in the 
masses of the third generation sleptons [13 ■ Since the effects are 
typically small, we leave a systematic investigation of extended 
scenarios for future work. 

We use the conservative BBN bounds considered in [2l| . The 
average EM energy release in one ri NLSP decay is assumed 
to be E-r/2, where E-r is the energy of the tau emitted in the 
dominant 2-body decay n Gr (cf. Fig. 16 of Ref. gl). With 
an EM energy release below Et/2, the light-shaded band can 
become smaller. For less conservative BBN constraints and/or 
enhanced EM energy release, the excluded mg region becomes 
larger. 



5 



CONCLUSION 

We provide the full SU(3)cXSU(2)lxU(1)y result for 
the relic density of thermally produced gravitino LSPs to 
leading order in the gauge couplings. Our result is ob- 
tained in a consistent gauge-invariant finite-temperature 
calculation and thus independent of arbitrary cutoffs. 
With this result, new gravitino and gaugino mass bounds 
emerge as a prediction of thermal leptogenesis. If super- 
symmetry is realized in Nature, these bounds will be ac- 
cessible at the LHC and the ILC. In particular, with a 
charged slepton NLSP, there will be the exciting possibil- 
ity to identify the gravitino as the LSP and to measure 
its mass. Confronting the measured gravitino mass with 
the predicted bounds will then allow for a unique test of 
the viability of thermal leptogenesis in the laboratory. 

We are grateful to T. Plehn, M. Pliimacher, G. Raffelt, 
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